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Figure  1.  Typical  SI  engine  envelope  of  gas 
temperature  and  pressure  histories  leading  up 
to  the  point  of  knock. 


[5].  Consequently,  existing  reduced  models  relevant  to 
autoignition  and  engine  knock  [6-10]  have  been  built  up 
from  skeleton  reaction  schemes  in  a  rather  empirical 
manner  to  include  representative  kinetic  steps,  which  can 
be  generic,  global  or  elementary  reactions.  Instead  of 
pursuing  the  details  of  oxidation  chemistry,  these 
models  have  been  developed  for  prediction  of  overall 
ignition  behavior.  For  paraffins  larger  than  C3,  a  typical 
autoignition  process  occurs  in  two  stages  characterized 
by  heat  release  as  in  Figure  2.  This  typical  two  stage 
ignition  characteristic  can  be  reproduced  by  the  reduced 
models.  Key  rate  parameters  are  generally  derived  by 
matching  experimental  results  (historically  ignitbn  delays) 
which  are  measured  in  idealized  facilities  such  as  rapid 
compression  machines  and  constant  volume  vessels. 

While  the  reduced  models  correctly  reproduced 
the  overall  trend  of  ignition  delays  as  a  function  of 
reaction  conditions  (temperature,  pressure,  and  reactant 
concentration)  and  fuel  structure,  significant  quantitative 
discrepancies  remain  between  the  model  prediction  and 
experimental  results.  Further,  the  empiricism  involved  in 
these  models  can  limit  the  scope  of  their  applications. 
To  enhance  their  predictive  capability,  reduced  models 
should  be  validated  against  more  detailed  information  on 
preignition  behavior  and  under  conditions  as  broad  and 
varied  as  possible. 

The  important  information  on  preignition  behavior 
that  has  been  overlooked  in  existing  models  includes 
the  preignition  heat  release  and  evolution  of  chemical 
species.  Experiments  have  shown  that  preignition 
reactions  can  significantly  increase  the  temperature  of 
end  gases  in  SI  engines  [3,111.  Further,  the  preignition 
heat  release  has  a  complex  dependence  witn  reaction 
conditions,  especially  the  temperature  [12,13].  Thus, 
accounting  for  preignition  heat  release  is  important  in 
reduced  models  in  order  to  enhance  their  predictive 
ability  under  broad  reaction  conditions.  Historically, 
reduced  models  have  not  been  validated  for  species 
prediction  and  do  not  perform  well  in  this  regard.  It  is 
highly  desirable  that  reduced  models  also  have  the 
ability  to  predict  at  least  key  chemical  species  or  species 
classes,  so  that  they  can  be  linked  and  tuned  with 
measured  species  data  to  better  represent  the  key 
chemical  processes 


Figure  2.  A  Epical  heat  release  profile  during  a  two 
stage  ignition  process. 


The  objective  of  the  current  study  is  to  develop  a 
reduced  chemical  kinetic  model  capable  of  predicting  the 
preignition  heat  release  and  key  chemical  species 
including  fuel  consumption,  formation  of  a  major 
oxidation  product,  CO,  and  production  of  other  species 
classes,  as  well  as  ignition  delays  at  low  and 
intermediate  temperatures.  This  work  is  a  continuatioii  of 
our  previous  efrort  in  which  we  extended  an  existing 
reduced  ignition  model  [9]  to  improve  its  prediction  of 
preignition  heat  release  arid  fuel  consumption,  and  built 
in  a  new  capacity  for  prediction  of  CO  formation  [14]. 
This  extended  model  was  tuned  to  match  our 
experimental  results  with  a  63  octane  blend  of  PRFs, 
63  PRF,  in  a  motored  engine  at  one  operating  conditbn. 
In  the  current  work,  we  modified  this  exterided  model 
using  experimental  data  from  our  recent  motored  engine 
study  of  n-heptane,  iso-octane  and  87  PRF  [15]. 
Modifications  were  made  to  improve  its  prediction  of 
preignition  heat  release,  fuel  consumption,  and  CO 
formatbn  and  to  account  for  productbn  of  other  species 
classes.  The  current  model  was  tuned  to  be  applic^le 
for  the  neat  PRFs,  87  PRF  and  63  PRF,  and  at  various 
engine  conditions. 

The  initial  model  which  we  extended  in  our  prior 
work  [14]  is  reviewed  next,  followed  by  a  detailed 
description  of  the  current  modified  model.  Results  of 
predbtions  using  the  current  model  are  then  presented 
and  discussed.  Rnally,  a  brief  discussbn  of  the  status  of 
this  model  is  provided. 


OVERVIEW  OF  PREVIOUS  MODEL 


In  our  previous  work  [14],  we  investigated  the 
ability  of  an  existing  model  of  18  reactions  and  13  active 
species  [9]  to  predict  heat  release  and  fuel  consumption 
during  the  preignition  reaction  of  63  PRF  in  a  motored 
research  engine.  The  baseline  model  is  listed  in  part  I 
of  Table  1 ,  where  RH  represents  the  fuel.  Given  the 
initial  fuel-air  mixture  concentration  and  temperature,  the 
chemical  kinetic  model  is  used  to  predict  temperature, 
heat  release  and  species  concentratbns  as  a  functbn  of 
time  or  crank  angb  by  integrating  the  coupled  rate  and 
energy  equations.  For  comparison,  we  independently 
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calculated  heat  release  from  measured  pressure  data 
using  a  standard  thermodynamic  model. 

We  found  that:  (i)  the  induction  time  of  heat 
release  can  be  matched  when  the  rate  parameters  of  the 
alkylperoxide  radical  (R02-)  isomerization  reaction  are 
chosen  between  the  values  suggested  for  iso-octane 
and  n-heptane;  (ii)  the  heat  release  predicted  by  the 
kinetic  model  is  about  20%  lower  than  that  calculated  from 
the  pressure  data;  and  (iii)  the  model  predicted  a  100% 
fuel  consumption  while  the  measured  fuel  consumption  is 
less  than  50%.  Analysis  of  the  baseline  model  showed 
that  in  its  present  form,  the  specific  heat  release  could 
not  be  matched. 

This  model  was  extended  by  adding  1 1  reactions 
and  7  active  species  to  account  for  the  oxidation  of 
aldehydes  (RCHO),  olefins  (C=C),  carbonyl  radicals 
(0=R-),  and  to  build  in  a  new  capacity  for  the  prediction 
of  CO  formation.  These  reactions  are  listed  in  Part  II  of 
Table  1.  The  heat  release,  fuel  consumption  and  CO 
formation  as  a  function  of  crank  angle  degree  (CAD) 
predicted  by  the  extended  model,  along  with  the 
predictions  by  the  baseline  model  and  experimental 
results,  are  presented  in  Rgure  3. 

The  results  show  that  the  extended  model 
significantly  improved  the  prediction  of  heat  release  and 
fuel  consumption.  In  addition,  the  extended  model 
predicted  CO  formation  reasonably  well.  However,  in 
spite  of  these  advances,  further  model  enhancements 
obviously  were  necessary  to  improve  these  predictions, 
and  to  allow  for  predictions  of  other  species  classes. 
Since  our  recent  experimental  study  [15]  provided 
detailed  species  information  for  n-heptane,  iso-octane 
and  87  PRF,  the  model  was  modified  and  adapted  for 
these  fuels.  The  final  version  of  the  model  was  also 
validated  against  the  experimental  results  of  63  PRF 
presented  in  Figure  3. 


RECENT  MODEL  MODIFICATIONS 


In  our  recent  motored  engine  study  of  autoignition 
chemistry  of  PRF’s,  we  mapped  the  overall  preignition 
reactivity  of  stoichiometric  n-heptane,  iso-octane,  and 
87  PRF  in  terms  of  exhaust  CO  production  as  a  function 
of  inlet  temperatures  (see  Figure  4).  Crank  angle 
resolved  evolution  profiles  of  fuel,  stable  intermediates 
and  oxidation  products  were  generated  for  lean  n- 
heptane  (0  =  0.2  and  0.3),  stoichiometric  iso-octane,  and 
stoichiometric  87  PRF  at  a  fixed  inlet  temperature  of 
376  K.  For  n-heptane  at  0  =  0.3,  the  measurements 
were  made  during  a  cyclically  repeatable  two  stage 
ignition  process  up  to  the  point  of  hot  ignition,  while  for 
n-heptane  at  0  =  0.2  and  other  fuels,  species  were 
measured  during  wcles  where  only  the  first  stage  of 
ignition  occurred.  The  results  have  been  reported  [15] 
except  for  n-heptane  at  0  =  0.2. 

The  experimental  results  of  n-heptane,  iso-octane, 
and  87  PRF  provided  a  basis  for  the  model  extension 
and  validation  in  the  current  study.  Besides  predicting 
heat  release,  fuel  consumption,  and  CO  formation,  the 
capability  to  predict  the  concentrations  of  other  species 
classes  was  also  a  target  of  this  work.  In  addition,  our 
overall  reactivity  mapping  results  at  various  engine 
operating  and  inlet  conditions  allowed  us  to  test  the 
dependence  of  the  model  prediction  on  reaction 
conditions  including  temperature  and  pressure  (or  charge 
density). 


Figure  3.  Comparison  of  the  heat  release,  fuel 
consumption  and  exhaust  CO  predictions  of  our 
previous  model,  along  with  the  predictions  by 
the  baseline  model  and  the  experimental  results 
with  63  PRF  [14]. 


Figure  4.  Exhuast  CO  molar  fractions  using 
stoichiometric  n-hetane,  iso-octane  and  87  PRF 
at  600  rpm  motored  engine  speed,  8.2 
compression  ratio,  and  108.2kPa  inlet  pressure 
[15]. 
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The  specific  model  modifications  were  made 
mainly  to  (i)  obtain  higher  specific  heat  release,  (ii)  allow 
prediction  of  other  species  classes,  and  (iii)  update 
conjugate  alkene  formation  submechanism.  Details  of 
the  model  development  follow. 

OBTAIN  HIGHER  SPECIFIC  HEAT  RELEASE  - 
As  Figure  3  shows,  at  the  end  of  the  first  stage  of 
ignition,  the  previous  model  predicted  fuel  consumption 
higher  than  the  experimental  result.  At  0  =  0.3  and 
376  K  inlet  temperature,  this  model  predicted  fuel 
consumption  too  high  in  comparison  with  experimental 
results  even  though  the  predicted  heat  release  profile 
matched  the  experimental  results.  Thus,  the  model 
needs  to  be  changed  to  obtain  higher  specific  heat 
release.  To  this  end,  the  oxidation  paths  of  carbonyl 
radicals  and  olefins  were  modified.  In  the  previous 
model,  the  oxidation  path  of  carbonyl  radicals  is 


C=R  +  C2  C=RC2-  (R1) 

which  is  reaction  28  in  Table  1 .  This  produces  a  high 
concentration  of  C=RC2-  at  the  end  of  reaction.  In 
reality,  these  large  radicals  oxidize  to  form  stable 
species  with  the  release  of  more  energy.  The  C=RC2- 
radicals  can  react  via  the  following  steps: 

C=RC2-  +  C=R02-  =>  2  C=RC-  +  C2  (R2) 

C=RC-  =>  RCHC  +  RCC-  (R3) 

if  the  fuel  is  iso-octane,  the  carbonyl  on  the  right  side  of 
R3  can  be  either  an  aldehyde  or  a  ketone.  An 
alternative  reaction  for  the  carbonyl  radical  is  the  direct 
decomposition  via  R4, 


0=R-  =>  Unsat.  Carbonyl  +  Rs-  (R4) 


where  Rs*  represents  alkyl  radicals  smaller  than  the 
parent  fuel.  For  iso-octane,  the  unsaturated  aldehydes 
are  methacrolein  and  acrolein.  When  the  fuel  is  n- 
heptane,  the  unsaturated  carbonyl  is  C=C=0.  Since 
C=C=0  has  never  been  experimentally  reported  in  the 
study  of  n-heptane,  this  path  is  probably  not  favored  for 
this  fuel. 

Since  the  reaction  paths  oxidizing  Rs*  and  RCO- 
to  stable  species  have  been  included  in  the  model,  the 
above  reactions  can  convert  0=R-  to  stable  species, 
resulting  in  higher  specific  heat  release.  However,  this 
would  introduce  three  additional  reactions  to  the  system. 
To  reduce  the  number  of  reactions  required,  the  result  of 
using  these  steps  was  compared  with  the  alternative  of 
using  a  single  global  step: 


0=R-  -I-  O2  =>  RCHO  -h  RCO-  1  O2  (R5) 
or 

0=R-  =>  RCHO  -I-  Rs-  (R6) 


No  significant  difference  in  model  predictions  was 
observed.  Without  affecting  the  overall  reaction  rate, 
either  R5  or  R6  can  even  be  combined  with  reaction  11 
in  Table  1,  OQ’O  =>  RCHO  +  0=R-,  into  one  global 
step: 


OQ'O  =>RCHO  -I-  RCHO  +  RCO-  -17.5  kcal/mol  (R7) 


or 

OQ’O  =>  RCHO  +  RCHO  +  Rs’  18.5  kcal/mol  (R8) 


It  should  be  noted  that  the  additional  O  atom  on  the 
right  side  of  R7  originates  from  R1.  As  mentioned 
above,  R4  is  not  favored  for  n-heptane,  and  therefore, 
R7  is  used  for  this  fuel.  While  both  R7  and  R8  are 
possible  for  iso-octane,  R8  gave  a  better  result  than 
K7.  Thus,  instead  of  pursuing  the  relative  importance  of 
these  two  pathways,  R8  was  selected  for  iso^octane  in 
the  current  study.  The  AH’s  of  R7  and  R8  are  derived 
based  on  the  structure  of  the  fuel  molecules  and  the 
chemistry  database  in  the  HCT  library  [16].  As  a  result 
of  these  simplifications,  the  number  of  total  reactions 
was  reduced  by  1 . 

Another  approach  to  obtain  higher  specific  heat 
release  is  to  modify  the  oxidation  paths  of  olefins.  In 
our  previous  model,  the  olefins  formed  by  Rs-  O2  <=> 
C=Cs  -I-  H02-  are  assumed  unreactive.  Model 
simulation  shows  that  the  concentration  of  C=Cs  is 
significantly  higher  than  that  of  oxygenates,  while  the 
experimental  results  showed  that  olefins  are  produced 
in  the  same  or  lower  concentration  than  oxygenates  for 
iso-octane  and  n-heptane,  respectively.  As  olefins  have 
higher  formation  enthalpy  than  oxygenates,  converting 
the  extra  C=Cs  into  oxygenates  can  result  in  additional 
heat  release.  The  oxidation  of  C=Cs  is  assumed  to 
follow  the  “Waddington”  mechanism  [17],  as  does  the 
olefin  C=C.  The  mechanism  proposes  that  after 
addition  of  OH-  to  the  C=C  double  bond,  the  resuitjng 
energetic  radical  -COH-C*-  will  be  stabilized  by  collision 
with  another  molecule.  Next,  O2  adds  to  the  stabilized 

radical  site  forming  -COH-COO-.  This  is  a  short-lived 
radical  which  decomposes  by  forming  a  six-membered 
transition  ring  and  splitting  into  two  carbonyls  and  OH-. 
However,  the  oxidation  paths  can  not  be  written  in  the 
same  form  as  for  C=C  (reactions  26  and  27  in  Table  1): 


C=Cs  -t-OH-  =>  -OHC-Cs-  26.5  kcal/mol  (R9) 
-OHC-Cs-  +  O2  =>2  RCHO  +OH-  -106  kcal/mol  (R1() 


The  reason  is  that  RCHO  will  be  oxidized  to  form 
C=Cs  via  reactions  19,  20  and  22  in  Table  1.  After  each 
cycle  the  number  of  C  atoms  is  reduced  by  1 .  When 
RCHO  becomes  a  formaldehyde,  C=Cs  can  not  be 
formed.  This  process  can  not  be  handled  by  a  reduced 
model  without  significantly  increasing  the  nutriber  of 
reactions.  Instead  of  pursuing  these  details,  the 
aldehydes  formed  through  this  path  are  labeled  as  Oxy 
(oxygenates)  and  assumed  unreactive.  Without  affecting 
the  model  performance,  these  two  reactions  can  also  be 
combined  into  one  global  reaction  as  R1 1 : 
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C=Cs  +  OH-  =>  2  Oxy  -f  OH- -76.5 kcal/mol  (R11) 


In  addition,  simulation  with  the  model  shows  that  the 
major  portion  of  the  olefins  are  C=Cs  formed  through 

the  smaller  alkyl  radicals  Rs’.  For  simplicity,  all  olefins 
are  treated  the  same  as  C=Cs.  The  subscript  “s”  can 
then  be  removed.  Thus,  the  number  of  reactions  is 
reduced  by  1.  The  rate  parameters  of  R11  are  selected 
from  the  work  of  Pitz  elal-  [18]  assuming  that  C=Cs  is 
iso-butene. 

ALLOW  PREDICTION  OF  OTHER  SPECIES 
CLASSES  -  There  are  several  considerations  in 
adapting  the  model  for  prediction  of  other  species 
classes.  First,  the  experimental  results  show  that,  a 
major  oxidation  product  class  for  PRF’s  is  the 
heterocyclic  ethers  which  can  account  for  30%  of  the  fuel 
reacted  for  iso-octane  and  26%  for  87  PRF.  Thus,  it  is 
desirable  to  include  this  reaction  path.  The  heterocyclic 
ethers  are  formed  through  two  steps: 


R02-  «  QOOH-7.  {R12) 

QOOH-7-=>  Ethers  +  OH-  (R13) 


ROOH  => 

RO-  + 

OH- 

43.6  kal/mol  (R16) 

RO-  => 

RCHO  + 

Rs- 

-10  kcal/mol  (R17) 

The  AH  and  rate  parameters  of  R16  are  the  same  as 
those  for  RsOOH  in  reaction  24  in  Table  1.  For  R17, 
the  rate  parameters  were  selected  based  on  the  work  of 
Westbrook  and  Pitz  [1 9],  while  the  AH  was  taken  from 
the  chemistry  database  in  the  HCT  library  [16]. 

Finally,  there  should  be  chemical  paths  for  all 
radicals  to  convert  to  molecular  species.  Thus,  reaction 
29  in  Table  1  was  removed  so  as  to  avoid  the  necessity 
for  additional  reactions  to  convert  RCOs-  into  stable 
products.  Such  detail  was  not  pursued  in  order  to 
minimize  the  total  number  of  reactions  in  the  model. 

After  these  modifications,  all  of  the  fuel  reacted  is 
eventually  converted  into  CO,  olefins,  heterocyclic 
ethers,  and  other  oxygenates,  which  is  consistent  with 
experimental  measurements. 

UPDATE  CONJUGATE  ALKENE  FORMATION 
SUBMECHANISM  -  In  the  previous  model,  conjugate 
alkenes  are  formed  via  direct  H  atom  abstraction  from 
alkyl  radicals  by  an  oxygen  molecule,  see  reactions  8  ard 
22  in  Table  1.  However,  recent  studies  [20,21]  have 
shown  that  conjugate  alkenes  are  primarily  formed 
through  coupled  mechanisms: 


the  -7  in  QOOH-7-  designates  that  the  ring  holds  seven 
atoms  including  oxygen  and  hydrogen  atoms.  For 
simplicity,  these  two  reactions  are  combined  into  a 
global  step: 

R02-  =>  Ethers  +  OH-  -25  kcal/mol  (R14) 

The  AH  of  R14  is  selected  based  on  the  chemistry 
database  in  the  HCT  library  [16].  For  the  rate  parameter 
selection,  the  activation  energy  of  R14  was  taken  from 
that  of  R12  in  a  detailed  model  [19]  assuming  the  parent 
fuel  is  iso-octane,  while  A-]  4  was  set  so  that  the  yield 
fraction  of  heterocyclic  ethers  match  the  experimental 
results. 

The  second  consideration  is  that  the  chain 
termination  reaction 

R-  +  R-  =>  RH  (R15) 

in  the  baseline  model  does  not  effectively  occur, 
because  neither  C14  nor  Cl  8  species  have  been 
reported  for  n-heptane  and  iso-octane.  This  reaction 
originated  from  the  model  developed  by  Halstead  el 
[6, 7]  in  which  R15  is  mathematically  necessary  to  limit  the 
reaction  rate  during  first  stage  ignition.  With  the 
extended  version  of  the  model,  this  reaction  is  not 
necessary  and  is  removed. 

Thirdly,  in  the  baseline  model,  ROOH  is  treated 
the  same  as  OQ’OOH.  For  species  prediction,  the  two 
species  need  to  be  differentiated,  because  OQ’OOH 
decomposes  to  eventually  form  a  RCHO,  an  0=R-  and 
an  OH-  while  ROOH  decomposes  to  produce  a 
RCHO,  a  Rs'  and  an  OH-  via  the  following  reactions: 


R"  +  O2 

<=> 

RO2* 

(R18) 

RO2* 

Olefin  -f  H02- 

(R19) 

RO2*  +  M 

R02-  +  M 

(R20) 

where  R02-*  and  M  represent  an  energized 
alkylperoxide  radical  and  any  third  body,  respectively. 
For  simplicity,  reactions  8  and  22  in  Tables  1  and  2  are 
replaced  by  the  following  two  reactions: 

R02’  =>  C=C 

+ 

H02-  17.5  kcal/mol 

(R21) 

Rs02‘  C=C 

+ 

H02-  4  kcal/mol 

(R22) 

The  activation  energy  of  R21  was  assigned  that  of  iso¬ 
octane  using  data  of  a  detailed  model  [19],  while  A22 
was  set  so  that  the  yield  fraction  of  conjugate  alkene 
matched  the  experimental  results.  It  was  successful  for 
iso-octane.  For  n-heptane,  however,  an  A  factor  of  R21 
which  matches  both  the  yield  fraction  of  conjugate 
alkenes  and  the  overall  reactivity  profile  could  not  be 
found.  For  simplicity,  the  A21  used  for  n-heptane  and 
87  PRF  was  taken  from  that  for  iso-octane.  This  resulted 
in  underprediction  of  conjugate  alkenes  for  these  two 
fuels  when  their  overall  reactivity  was  matched.  The  AH 
of  R21  was  taken  from  the  chemistry  database  in  the 
HCT  library  [16].  The  rate  parameters  of  R22  were 
taken  from  a  detailed  mechanism  of  ethane  used  by 
Wagner  et  a].  [22],  while  the  AH  for  this  reaction  is 
derived  from  reactions  21  and  22  in  Table  1.  The  new 
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conjugate  alkene  formation  mechanism  improved  the 
sensitivity  of  the  model  to  charge  density. 

One  additional  modification  was  the  addition  of 
R23,  which  is  believed  important  for  alkanes  containing 
fewer  than  4  carbon  atoms.  The  AH  and  rate  parameters 
were  taken  from  reaction  1 7  in  Table  1 . 

RH  +  Rs02'  =>  RsOOH  +  R-  8 kcal/mol  (R23) 

The  new  version  of  our  model  is  detailed  in  Table 
2.  The  total  number  of  reactions  and  number  of  active 
species  are  the  same  as  in  Li  Mai-  [14],  but  the  details 
are  different.  Results  using  this  model  are  presented 
next.  In  the  following  section,  the  comparison  between 
the  new  model  and  experiments  is  presented,  along 
with  a  discussion  of  the  rate  parameters  recommended 
for  key  reactions. 


RESULTS  WITH  THE  CURRENT  MODEL 


Two  sets  of  experimental  data  were  used  to  test 
and  tune  the  model.  The  first  set  is  the  crank  angle 
resolved  heat  release  and  species  data  at  376  K  inlet 
temperature  for  lean  n-heptane  (o  =  0.3  and  0.2)  and  for 
stoichiometric  iso-octane  and  87  PRF.  The  second  set  is 
the  exhaust  CO  mapping  data  at  various  inlet 
temperatures  and  compression  ratios  for  iso-octane  and 
87  PRF. 

The  modeling  procedures  were  described  in  our 
revious  paper  [14],  but  are  briefly  reviewed  here, 
irst,  core  temperatures  and  heat  release  are  calculated 
from  pressure  data  using  a  thermodynamic  model. 
Secondly,  the  energy  equation  is  integrated  to  simulate 
compression  temperature  assuming  no  chemical 
reactivity.  The  heat  transfer  coefficient  is  set  so  that  the 
simulated  temperatures  match  the  core  temperature 
profiles  calculated  from  pressure  data  before  sensible 
preignition  heat  release  is  observed.  Figure  5  shows 
the  example  match  between  the  pressure  derived  and 
model  simulated  temperature  profiles  at  376  K  inlet 
temperature  where  detailed  species  data  were 
generated  for  n-heptane,  iso-octane  and  87  PRF. 


Figure  5.  Core  temperature  profiles  derived  from 
pressure  data  and  the  temperature  profiles 
simulated  by  the  model  assuming  no  reactivity. 


After  setting  the  heat  transfer  coefficient,  the  chemical 
kinetic  rate  equations  and  the  energy  equation  are 
integrated  to  predict  temperature  and  species 
concentrations  as  a  function  of  crank  angle  degree. 

In  this  study,  an  effort  was  made  to  account  for  the 
effects  of  residual  gases  on  model  prediction.  Both 
experiments  and  detailed  chemistry  modeling  have 
shown  that  the  partial  oxidation  products  in  the  residual 
gases  can  promote  reactivity  [23,24].  Westbrook  et  §1- 
[25]  tested  the  effects  of  pro-knock  additives  on  the 
ignition  time  using  a  detailed  chemical  kinetic  rr»del.  The 
tested  additives  include  alkyl  hydroperoxides,  hydrogen 
peroxide,  aldehydes,  olefins,  etc.  It  was  found  that  the 
alkyl  hydroperoxides  are  the  most  effective  in  reducing 
the  ignition  time  because  they  can  decompose  to  form 
the  active  radical,  OH-.  In  the  reduced  model,  the  most 
important  chain  branching  agents  are  OQ’OOH.  Thus, 
the  initial  concentration  of  this  species  was  taken  into 
account.  The  procedures  are  as  follows.  First,  assuming 
the  initial  charge  solely  consists  of  fuel  and  air, 
representing  the  first  motored  cycle,  the  concentrations 
of  these  species  are  calculated  using  the  reduced 
model.  The  calculated  concentration  of  OQ’OOH  and 
the  residual  fraction  03%)  were  used  to  determine  the 
initial  concentration  of  these  species  in  the  next  motored 
cycle  simulation.  Such  procMures  were  repeated  until 
the  initial  concentration  of  OQ’OOH  remained 
unchanged. 

MODELING  CAD  RESOLVED  DATA  -  The 
measured  CAD  resolved  data  of  iso-octane,  n-heptane 
and  87  PRF  at  376  K  inlet  temperature  were  cornpared 
with  model  predictions.  The  simulation  started  with  Iso¬ 
octane.  Using  the  parameters  suggested  in  the 
baseline  model,  the  reactivity  was  calculated  to  start  later 
in  the  cycle  than  measured  experimentally.  Thus,  rate 
parameters  needed  to  be  calibrated  to  match  the 
observed  induction  time.  In  simulating  knock  in  an  engine 
using  the  Hu  and  Keck  model  [9],  Cowart  Mai-  126] 
calibrated  the  model  by  reducing  the  activation  energy  of 
the  R02-  isomerization  reaction  (reaction  3  in  Tables  1 
and  2): 


R02-  QOOH-  (R24) 


from  22.4  to  21 .4  kcal/mol.  In  addition,  recent  ignition 
delay  measurements  by  Park  M  ll-  [27]  are  significantly 
shorter  than  the  measurements  in  the  Thornton  rapid 
compression  machine  to  which  the  Hu  and  Keck  model 
was  calibrated  [9].  Thus,  the  rate  parameters  would 
need  to  be  adjusted  to  make  the  reaction  rate  faster, 
which  is  consistent  with  this  study.  The  parameter 
adjusted  here  is  the  activation  energy  of  reaction  7,  which 
was  reduced  from  43  to  40  kcal/mol  to  match  the  induction 
time  of  iso-octane.  This  parameter  was  selected 
because  there  is  relatively  high  uncertainty  in  the  original 
value.  Another  important  reason  for  selecting  this 
reaction  is  that  it  affects  ignition  delay  while  having  only  a 
small  effect  on  the  magnitude  of  the  preignition  heat 
release.  The  same  activation  energy  for  this  reaction  was 
used  for  both  87  PRF  and  n-heptane. 

The  most  effective  reaction  in  modifying  the  heat 
release  profile  is  R25  (reaction  17  in  Table  2): 
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RCHO  +  OH-  =>RCO-  +  H2O 


(R25) 


CAD 


Figure  6.  Comparison  of  heat  release,  fuel  consumption 
and  [CO]  predictions  of  the  current  model  with 
the  experimental  results  of  iso-octane  at  8.2 
compression  ratio,  600  rpm,  108.2  kPa  inlet 
pressure,  and  376  K  inlet  temperature. 


This  reaction  has  little  impact  on  the  induction  time  of  the 
first  stage  of  ignition.  The  baseline  rate  parameters  of 
this  reaction  were  taken  from  the  HCT  database  [16] 
assuming  that  RCHO  is  acetaldehyde.  To  match  the 
heat  release  profile.  A-)  g  had  to  be  set  at  3.7  times  the 
baseline  value.  The  rationale  for  making  this  adjustment 
is  explained  shortly  with  the  results  of  modeling  n- 
heptane  oxidation.  For  other  reactions,  A-|4  and  A21 
were  set  to  match  the  yield  fraction  of  heterocyclic  ethers 
and  conjugate  alkenes,  which  are  30%  and  9%, 
respectively. 

The  results  of  heat  release,  fuel  consumption,  and 
CO  concentration  predicted  by  the  model,  along  with 
the  experimental  data,  are  shown  in  Figure  6.  It  can  be 
seen  that  the  computed  profiles  are  quite  close  to  the 
experimental  results.  For  other  species  classes,  the 
predicted  concentration  of  oxygenates  other  than 
heterocyclic  ethers  is  approximately  the  same  as  olefins 


at  the  end  of  reaction,  which  is  consistent  with 
experiments. 

The  simulated  and  experimental  results  of  n- 
heptane  at  0  =  0.3  and  0.2  are  presented  in  Figure  7. 
As  suggested  by  Hu  and  Keck  [9],  the  rate  parameters 
of  R02-  isomerization  (R24)  are  the  major  variables  to 
account  for  the  fuel  structure  effects  on  ignition  delay. 
Using  the  same  values  of  A24  and  E24  in  the  baseline 
model,  the  induction  times  were  matched  for  both 
equivalence  ratios.  However,  the  magnitude  of  first 
stage  heat  release  was  lower  than  in  the  experiments. 
In  order  to  match  the  preignition  heat  release,  the  A25 
was  tuned  to  a  value  1.65  times  the  one  in  the  HCT 
database  [16]  assuming  that  RCHO  is  acetaldehyde. 
The  third  rate  parameter  adjusted  was  A14  to  match  the 
yield  fraction  of  heterocyclic  ethers.  The  value  selected 
IS  3.3  times  that  of  iso-octane.  With  these  rate 
parameter  adjustments,  the  predicted  profiles  of  heat 
release,  fuel  consumption  and  CO  formation  are  close  to 
experimental  results  except  that  the  second  stage 
ignition  delay  occurred  4  CAD  earlier  than  in  the 
experiment  with  0  =  0.3.  Consequently,  the  A  factor  of 
the  chain  branching  reactbn: 


H2O2  +  M  =>  2  OH-  +  M  (R26) 


(reaction  10  in  Table  2)  was  multiplied  by  0.6  to  obtain  a 
better  match  of  the  ignition  delay  of  the  second  stage 
ignition.  The  predicted  and  experimental  results  are 
shown  in  Figure  7.  For  other  species  classes,  the 
predicted  concentration  of  oxygenates  other  than 
heterocyclic  ethers  is  about  30%  higher  than  olefins  at  the 
end  of  the  first  stage  but  before  the  second  stage 
ignition,  which  is  again  consistent  with  experiments. 

The  modeling  results  of  n-heptane  are  exciting  in 
that,  with  a  value  of  A25  close  to  that  of  the 
corresponding  elementary  reaction,  the  heat  release 
profile,  fuel  consumption,  and  CO  formation  at  both 
equivalence  ratios  are  in  good  agreement  with 
experiments.  However,  it  is  necessary  to  justify  the  use 
of  a  larger  A25  for  iso-octane  than  n-neptane.  In 

addition  to  leading  to  CO  formation,  R25  removes  OH- 
and  therefore  reduces  preignition  reactivity.  In  reality, 
olefins  may  also  play  such  a  role.  This  is  particularly  true 
for  iso-octane  because  its  major  intermediate  species, 
iso-butene,  reduces  reactivity  at  low  temperatures. 
However,  for  simplicity,  the  current  model  assumes  that 
all  olefins  follow  only  the  “Waddington"  mechanism, 
resulting  in  an  exclusion  of  this  possible  inhibition  role. 
It  is  expected  that  A25  is  similar  for  both  n-heptane  and 
iso-octane  if  the  inhibition  path  of  iso-butene  is  taken 
into  account.  But  this  would  significantly  increase  the 
number  of  total  reactions.  Hence,  these  details  were 
not  taken  into  account  in  the  current  study. 

For  87  PRF,  only  those  rate  parameters  which 
are  different  for  n-heptane  and  iso-octane  were  tuned. 
For  R24,  the  activation  energy  for  iso-octane  was  used 
for  87  PRF  and  A24  was  set  4  times  that  of  iso-octane  to 
match  the  observed  induction  time.  The  A25  selected  to 
match  the  heat  release  profile  was  3.5  times  the  one  in 
the  HCT  database  [16]  assuming  that  RCHO  is 
acetaldehyde.  In  order  to  match  ethers  yield  fraction 
(26%),  the  A-|4  used  was  2.5  times  that  for  iso-octane. 
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Figure  7.  Comparison  of  heat  release,  fuel  consumption 
and  [CO]  predictions  of  the  current  model  with 
the  experimental  results  of  n-heptane  at  8.2 
compression  ratio,  600rpm,  108.2  kPa  inlet 
pressure,  and  376  K  inlet  temperature. 


The  results  of  heat  release,  fuel  consumption,  and  CO 
concentration  predicted  by  the  model,  along  with  the 
experimental  data,  are  shown  in  Figure  8.  The  heat 
release  profile  is  matched  very  well,  but,  the  CO  and 
fuel  consumption  differ  from  the  experiments  by 
approximately  20%.  For  other  species  classes,  the 
predicted  concentration  of  oxygenates  other  than 
heterocyclic  ethers  is  about  20%  higher  than  olefins  at  the 
end  of  reaction,  while  the  experimental  results  show  that 
the  former  is  about  80%  higher  than  the  latter.  Thus,  the 
species  prediction  for  87  PRF  is  not  as  good  as  for  the 
neat  PRF’s. 

The  above  predictions  were  conducted  at  one  inlet 
condition.  The  performance  of  the  model  at  other 
engine  operation  conditions  is  examined  next. 

MODELING  AT  VARIOUS  ENGINE 
CONDITIONS  --  The  major  physical  factors  affecting 
oxidation  reactivity  at  low  and  intermediate  temperatures 
are  temperature  and  pressure  (or  charge  density).  It  is 


Figure  8.  Comparison  of  heat  release,  fuel  consumption 
and  [CO]  predictions  of  the  current  model  with 
the  experimental  results  of  87  PRF  at  8.2 
compression  ratio,  600  rpm,  108.2  kPa  Inlet 
pressure,  and  376  K  inlet  temperature. 


important  that  the  model  has  the  ability  to  account  for  the 
effects  of  these  factors. 

As  shown  in  Figure  4,  as  the  inlet  temperature 
was  increased,  the  reactivity  of  iso-octane  first  increased, 
followed  by  a  plateau,  and  then  decreased,  exhibiting  a 
negative  temperature  coefficient  (NTC)  behavior.  No 
autoignition  occurred.  The  reactivity  of  87  PRF  was  much 
higher  than  iso-octane,  but  its  profile  as  a  function  of 
temperature  is  similar  to  iso-octane,  ^oradic 
autoignitbn  (two  stage  ignition)  was  observed  as  the 
inlet  temperature  was  increased  to  377  K.  However, 
when  the  inlet  temperature  reached  480  K,  autoignition 
ceased.  As  the  inlet  tenrperature  was  further  increased 
to  500  K,  autoignition  still  did  not  occur  and  exhaust  CO 
was  further  decreased.  This  observation  suggests  that 
the  reactivity  or  magnitude  of  preignition  reaction  has  a 
significant  effect  on  the  second  stage  of  ignition. 

The  existing  general  hydrocarbon  oxidation 
mechanism  can  explain  the  NTC  behavior  in  terms  of  a 
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shift  of  R02-  reaction  mechanism  to  alkene  formation 
paths.  This  mechanism  is  included  both  in  the  baseline 
model  and  in  the  extended  reactions.  Thus,  the  model 
can  at  least  qualitatively  account  for  the  effects  of 
temperature.  The  ability  of  the  model  to  quantitatively 
predict  the  NTC  behavior  was  tested  using  the  overall 
reactivity  mapping  results  presented  in  Figure  4. 

As  to  the  pressure  effects,  experiments  show  that 
preignition  reactivity  increases  as  the  pressure  (or  charge 
density)  increases.  Benson  [28]  qualitatively  analyzed 
the  effects  of  pressure  on  the  temperature  at  which  the 
R02-  shift  from  degenerate  chain  branching  paths  to 
conjugate  alkene  formation  paths.  The  analysis  is  based 
on  the  following  reactions: 


R-  + 

02 

o  R02' 

{R27) 

R02'  + 

RH 

=>  ROOH  + 

R' 

(R28) 

R-  + 

02 

o 

II 

o 

+ 

H02' 

{R29) 

Based  on  the  steady-state  analysis  and  assuming  that 
R02-  and  R-  are  in  equilibrium  via  R27,  the  rates  of 
production  of  ROOH  and  the  olefin  is  developed  and 
found  to  give  the  ratio 


d[C— C]  _  k29 

cl[ROOH]  "  k28k27[RH] 


Thus,  at  the  same  temperature,  when  the  charge  density 
is  increased,  the  concentration  of  the  fud  [RH]  is 
increased  and  the  formation  of  ROOH  is  favored, 
resulting  in  higher  reactivity.  Since  ROOH  is  a  major 
chain  branching  agent  for  light  HC’s,  equation  (1)  applies 
for  light  HC’s  only.  For  larger  HC’s  like  PRFs,  the  major 
chain  branching  paths  are  through  R24.  The  charge 
density  does  not  directly  affect  the  competition  between 
chain  branching  agent  formation  and  olefin  formation. 
Thus,  the  baseline  model  has  a  relatively  weak 
dependence  of  reactivity  on  charge  density.  For  Rs- 
reactions  in  the  previous  model,  however,  Benson’s 
analysis  applies.  With  the  updated  conjugate  formation 
path  R22  in  the  current  model,  the  rate  of  production  of 
RsOOH  and  the  olefin  is 


cl[C— Cs]  _  k22 

d[RsOOH]  "  k23  [RH] 


It  should  be  noted  equation  (2)  holds  without  the 
assumption  of  the  equilibrium  between  R02-  and  R-. 
Thus,  the  extended  reactions  in  the  current  model 
improve  the  pressure  dependence  of  the  model. 

To  test  the  model’s  sensitivity  to  pressure,  or 
charge  density,  an  exhaust  CO  mapping  of  87  PRF  at 
CR  =  5.2  was  also  compared  with  model  predictions.  In 
that  experiment,  the  engine  was  operated  at  900  rpm, 
165  kPa  with  inlet  pressure  and  the  charge  density  at  top 
dead  center  (TDC)  being  about  15%  lower  that  at  CR  = 
8.2. 


Figure  9.  Comparison  of  the  exhaust  CO  predictions  of 
the  current  model  with  the  experimental  results 
of  87  PRF  and  iso-octane  at  various  engine 
conditions. 


The  experimental  results  and  model  predictions 
are  presented  in  Figure  9.  It  can  be  seen  that  the  NTC 
behavior  of  iso-octane  and  87  PRF  were  reproduced. 
The  predicted  dependence  of  reactivity  on  pressure,  or 
charge  density  of  87  PRF  was  also  in  agreement  with 
experiments. 

During  the  reactivity  test  of  87  PRF,  sporadic 
autoignition  was  observed  experimentally  at  inlet 
temperatures  between  377  and  480  K,  while  the  rnodel 
predicted  no  autoignition  using  the  averaged  iriitial 
parameters.  It  is  desirable  that  the  sporadic  autoignition 
phenomena  be  interpreted  by  the  reduced  model.  In 
this  regard,  the  causes  for  this  phenomena  were 
analyzed  in  order  to  identify  the  initial  parameters  which 
need  to  be  changed  for  modeling  the  autoignition  cycle. 
From  the  pressure  traces,  it  was  observed  that  there  is 
no  apparent  reactivity,  or  heat  release  in  the  cycle 
before  an  autoignition  cycle.  The  model  simulation 
shows  that  for  a  cycle  at  the  transition  state  from  minimal 
to  significant  reactivity,  the  level  of  OQ’OOH  can  be  up 
to  hundreds  of  ppm.  On  the  other  hand,  for  cycles  with 
sig^nificant  reactivity,  the  simulated  concentration  of 
OQ’OOH  is  only  about  0.5  ppm.  Thus,  it  is  highly  likely 
that  the  “non-reactivity’'  cycle  provided  a  high 
concentration  of  active  species  which  promoted  reactivity 
and  produced  autoignition  in  the  next  cycle.  This 
hypothesis  is  supported  by  the  experimental  work  of 
Walcutt  filai-  [29],  in  which  they  found  that  the 
thermoneutral  reactions  occurring  prior  to  the  onset  of 
cool  flames  result  in  the  decrease  in  the  anti-knock 
performance  of  the  fuel. 

Following  the  above  analysis,  the  initial  level  of 
OQ’OOH  in  the  model  was  gradually  increased  at  the 
inlet  temperature  of  377  K,  where  autoigriition  was  first 
experimentally  observed.  The  simulation  predicted 
autoignition  when  the  initial  level  of  OQ’OOH  reached 
130  ppm.  This  indicated  that  if  the  “non-reactivity"  cycle 
produces  1000  ppm  OQ’OOH,  which  is  possible,  the 
next  cycle  will  autoignite.  Using  the  same  level 
OQ’OOH,  the  model  did  not  predict  autoignition  for 
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Figure  10.  Comparison  of  the  heat  release,  fuel 
consumptions  and  CO  predictions  of  the  current 
model  with  the  experimental  data  of  63  PRF  at 
5.2  compression  ratio,  900  rpm,  165  kPa  inlet 
pressure,  and  445  K  inlet  temperature. 


87  PRF  at  lower  temperatures  nor  autoignition  of  iso¬ 
octane  at  any  temperature,  which  is  consistent  with 
experiments.  Thus,  this  sporadic  autoignition  of  87  PRF 
can  be  interpreted  by  the  reduced  model.  However, 
due  to  the  uncertainties  in  the  initial  level  of  OQ’OOH,  a 
quantitative  prediction  on  the  inlet  temperature  range 
where  autoignition  occurred  is  impossible. 

Since  significant  modifications  have  been  made  to 
our  previous  model,  it  is  prudent  to  test  the  current 
model  against  the  experimental  results  of  63  PRF 
presented  in  Figure  3.  Following  the  same  method  for 
modeling  87  PRF,  A24  was  set  to  match  induction  time, 
while  A25_  was  set  to  match  the  magnitude  of  heat 
release.  Since  the  hetrocyclic  ethers  were  not  reported 
for  that  fuel,  the  value  of  A14  was  assigned  by  a  RON 
based  linear  interpolation  using  the  Ai4’s  of  n-heptane 
and  87  PRF.  The  final  value  of  A24  is  9.5  times  that  of 
iso-octane,  while  A25  is  2.7  times  that  in  the  HCT 
database  [16]  assuming  that  RCHO  is  acetaldehyde. 
The  results  of  heat  release,  fuel  consumption,  and  CO 
concentration  predicted  by  the  model,  along  with  the 
experimental  data,  are  shown  in  Figure  10.  It  can  be 


seen  that  the  predicted  heat  release  profile  is  much 
improved  compared  with  Figure  3.  The  fuel 
consumption  and  CO  predictions  differ  from 
experiments  by  about  22  and  25%,  respectively.  It  is 
also  observed  that  the  problem  with  continued  CO 
formation  at  late  CAD’s  predicted  by  the  previous 
model  disappears. 


DISCUSSION 


Table  2  represents  the  final  reduced  mechanism 
and  recommended  rate  parameters,  including  the 
modifications  necessary  to  accommodate  different  fuels 
structure.  As  expected,  the  fuel  specific  parameters 
change  monotonically  with  the  octane  number  of  the 
fuels.  In  addition,  using  the  rate  parameters  developed 
based  on  the  CAD  resolved  measurements  at  one  inlet 
engine  condition,  the  model  reproduced  the  complicated 
oxidation  behavior  at  various  engine  inlet  and  operating 
conditions.  This  behavior  indicates  that  the  modification 
of  these  rate  parameters  was  not  just  an  empirical 
exercise,  but  one  that  captured  the  essence  of  the 
underlying  chemical  processes. 

In  comparison  with  other  reduced  models,  two 
improvements  were  made  with  the  current  model: 

(1)  The  mechanisms  (primarily  via  secondary 
.  reactions)  necessary  to  account  for  preignition 

heat  release  have  been  included. 

(2)  All  of  the  fuel  reacted  is  eventually  converted 
to  molecular  oxidation  products  which  can  be 
divided  into  CO,  olefins,  heterocyclic  ethers 
and  other  oxygenates. 

Under  the  experimental  conditions  examined,  the 
model  predictions  agreed  fairly  well  with  experiments. 
Nonetheless,  there  are  a  number  of  factors  that  must  be 
considered  prior  to  full  utilization  of  this  model  for  engine 
applications.  First,  the  model  needs  to  be  tested  under 
broader  engine  conditions,  particularly  at  higher 
pressures  and  higher  engine  speed,  so  that  the  fuel/air 
mixture  experiences  conditions  closer  to  those  of  the 
end  gas  in  fired  engines. 

Second,  the  number  of  reactions  needs  to  be 
optimized.  Existing  reduced  models  consist  of  reactions 
ranging  from  a  few  to  nearly  a  hundred.  From  the 
viewpoint  of  engine  applications,  a  model  with  the 
sufficient  prediction  ability  and  the  smallest  number  of 
reactions  is  desirable.  The  issue  regarding  this  trade-off 
remains  unaddressed,  and  definitely  deserves  additional 
effort.  We  feel  that  there  is  a  potential  for  reducing  the 
number  of  reactions  in  the  current  model  without 
significantly  affecting  its  predictions. 

Third,  the  model  needs  to  be  adapted  for  a 
broader  range  of  fuels.  Since  most  paraffins  follow 
similar  chemical  pathways,  the  present  model 
developed  for  PRF’s  can  be  conveniently  adapted  to 
other  paraffins.  However,  in  order  to  develop  a  model 
for  practical  fuels,  the  model  needs  to  be  extended  to 
include  initiation  steps  for  olefins,  aromatics,  oxygenates 
and  their  mixtures. 
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SUMMARY  AND  CONCLUSIONS 


We  developed  a  reduced  chemical  kinetic  model 
for  the  prediction  of  major  oxidation  behavior  of  PRF’s 
at  low  and  intermediate  temperatures,  including  ignition 
delay,  preignition  heat  release,  fuel  consumption,  CO 
formation  and  production  of  other  species  classes.  The 
model  development  was  carried  out  in  two  phases.  In 
the  first  phase,  which  was  previously  reported  [14],  an 
existing  reduced  ignition  model  of  18  reactions  and 
1 3  active  species  was  extended  to  account  for  specific 
heat  release  and  CO  formation,  ending  up  with  7  new 
active  species  and  1 1  additional  reactions.  This  model 
was  tuned  to  match  the  experimental  results  of  a 
63  PRF,  at  one  engine  operating  condition. 

In  the  second  phase,  the  model  was  modified 
utilizing  n-heptane,  iso-octane  and  87  PRF  experimental 
data  from  our  recent  motored  engine  study  of  PRF 
autoignition  chemistry  [15],  Modifications  primarily 
includfed  the  refinement  of  carbonyl  and  olefin  reactions 
for  higher  specific  heat  release,  incorporation  of  reaction 
paths  to  generate  key  classes  of  species,  and  updating 
the  conjugate  alkene  formation  submechanism.  The  total 
number  of  active  species  and  reactions  remains  the 
same  as  the  initially  extended  model,  but  with  different 
details.  Accommodating  the  different  fuels  requires 
mainly  adjusting  the  rate  parameters  of  the  R02- 
isomerization  reaction,  the  reaction  of  aldehydes  with 
OH-,  and  the  reaction  of  forming  cyclic  ethers  for  each 
fuel. 

At  the  motored  engine  condition  where  detailed 
species  data  were  generated,  the  model  can  reproduce 
the  ignition  delay  and  the  magnitude  of  the  preignition 
heat  release.  Fuel  consumption  and  CO  formation 
predictions  differed  from  experiments  by  at  most  25% 
for  all  of  the  four  fuels.  Predictions  for  other  species 
classes  generally  agreed  with  experiments.  The  overall 
reactivity  measurements  at  various  inlet  conditions  were 
also  used  to  test  the  dependence  of  the  model 
predictions  on  temperature  and  charge  density.  The 
NTC  behavior  of  iso-octane  and  87  PRF  was 
reproduced  by  the  model.  In  addition,  the  lower 
reactivity  of  87  PRF  at  a  lower  compression  ratio  was 
also  reproduced,  indicating  that  the  model  can  account 
for  the  effects  of  pressure  or  charge  density. 

These  results  indicated  that  the  current  reduced 
kinetic  model  is  a  useful  tool  for  prediction  of  major 
oxidation  behavior  at  low  and  intermediate 
temperatures. 
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Tablel 

Previous  Reduced  Chemical  Kinetic  Model 
Part  I:  The  Baseline  Model  [9] 


A.  13  Active  Species 


1. 

RH 

2. 

02 

3. 

R- 

4. 

R02- 

5. 

QOOH 

6. 

OOQOOH- 

7. 

OQ’O- 

8. 

OH- 

9. 

H02- 

10. 

HOOH 

11. 

OQ’OOH 

12. 

RCHO 

13. 

C=C 

B.  1 8  Reactions  (  unit:  cc,  mole,  s,  kcal ) 

Arrhenius  parameters  of  rate  constants  k  =  A  e  "E/RT  fQ,-  heptane  at  600  K<T<1100K 


Equilibrium  k+  k" 


Reaction 

AH®3(X) 

log  A 

E 

log  A"*" 

E  + 

log  A’ 

E- 

1 

RH  +  O2  =>  R-  -H  H02- 

46.4 

1.5 

46.0 

13.5 

46.0 

12.0 

0.0 

2 

R-  +  O2  <=>  R02' 

-31.0 

-1.4 

-27.4 

12.0 

0.0 

13.4 

27.0 

3 

R02-  <=>  QOOH-  (n-heptane) 

7.5 

0.9 

8.0 

11.9 

19.0 

11.0 

11.0 

(iso-octane) 

7.5 

0.0 

11.4 

11.0 

22.4 

11.0 

11.0 

4 

QOOH-  +  O2  «  OOQOOH- 

-31.0 

-1.9 

-27.4 

11.5 

0.0 

13.4 

27.4 

5 

OOQOOH-  =>  OQ’OOH  +  OH- 

-26.6 

11.3 

17.0 

6 

OH-  -f  RH  =>  H2O  +  R- 

-23.5 

13.3 

3.0 

7 

OQ’OOH  =>  OQ’O-  +  OH- 

43.6 

15.6 

43.0 

8 

R-  +  O2  <=>  C=C  +  HO2- 

-13.5 

0.0 

-13.5 

11.5 

6.0 

11.5 

19.5 

9 

HO2-  -f  HO2-  =>  HOOH  +  O2 

-38.5 

12.3 

0.0 

10 

HOOH  +  M  =i>  20H-  -I-  M 

51.4 

17.1 

46.0 

11 

OQ’O-  =>  RCHO  +  0=R- 

8.5 

14.0 

15.0 

12 

QOOH-  =>C=C  +  RCHO  +  OH- 

-3.0 

14.4 

31.0 

13 

RO2-  RCHO  =>  ROOH  +  RCO- 

-0.6 

11.45 

8.6 

14 

HO2-.+  RCHO  =?.  HOOH  +  RCO- 

-0.6 

11.7 

8.64 

15 

C=C  +  HO2-  =>  Epox  +  OH- 

-0.23 

10.95 

10.0 

16 

HO2-  4  RH  «  R-  4  HOOH 

8.0 

0.9 

8.0 

11.7 

16.0 

10.8 

8.0 

17 

R02-  4  RH  »  ROOH  4  R- 

8.0 

1.1 

8.0 

11.2 

16.0 

10.1 

8.0 

18 

R-  4  R-  RH 

-85.0 

13.2 

0.0 

Note:  ROOH  and  OQ’OOH  are  treated  as  identical;  0=R-,  RCO-  and  Epox  are  assumed  inactive. 
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Tablel  (cont.) 

Previous  Reduced  Chemical  Kinetic  Model 
Part  II:  Extended  Reactions  [14] 


A.  7  Activ 


1.  RCO-  2.  Rg-  3.  Rs02-  4.  RgOOH  5.  RgO-  6.  -OHC-C-  7.  RCOa* 


B.  1 1  Reactions  (  unit:  cc,  mole,  s,  kcal ) 

Arrhenius  parameters  of  rate  constants  k  =  A  e  'E/RT 

Equilibrium 

Reaction  AH^soo  log  A  E 

19  RCHO  +  OH-  =>  RCO-  +  H2O  -33.5 

20  RCO-  +  M  Rg-  +  CO  +  M  8.7 

21  Rg- +  O2  <=>  Rs02-  -31.0  -1.4  27 j 

22  Rg-+  O2  <=>  C=Cg  +  H02-  -13.5  0.0  13.1 

23  RCHO  +  Rg02-  =>  RgOOH  -f  RCO-0.6 


RgOOH 


RgO-  +  OH- 


RgO-  +  O2  =>  Rs'O  +  H02- 
C=C  +  OH-  =j.  -OHC-C  - 
-OHC-C-  +  O2  =>  2  RCHO 
0=R-  +  O2  <=>  0=R02- 
RCO-  +  O2  <=>  RCO3- 


AH°3(x) 

log  A 

E 

log  A+ 

E  + 

log  A' 

E- 

-33.5 

12.62 

0.52 

8.7 

16.78 

15.0 

-31.0 

-1.4 

27.4 

12.0 

0.0 

13.4 

27.0 

-13.5 

0.0 

13.5 

11.5 

6.0 

11.5 

19.5 

)-0.6 

11.45 

8.6 

43.6 

15.6 

43.0 

-26.5 

10.6 

2.14 

26.5 

12.68 

-0.79 

■106.0 

10.0 

-4.0 

-31.0 

-1.9 

-27.4 

11.5 

0.0 

13.4 

27.4 

-31.0 

-1.9 

-27.4 

11.5 

0.0 

13.4 

27.4 

Note:  Rg'O  and  C=Cg  are  assumed  inactive. 
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Table  2 

Current  Reduced  Chemical  Kinetic  Model  --  Part  1 


A.  13  Active  Species 

I.  RH  2.  O2  3.  R-  4.  RO2’  5.  QOOH- 

6.  OOQOOH-  7.  OQO-  8.  OH-  9.  H02-  10.  HOOH 

II.  OQ'OOH  12.  RCHO  13.  C=C 

B.  16  Reactions  ( unit:  cc,  mole,  s,  kcal ) 

Arrhenius  parameters  of  rate  constants  k  =  A  e  "E/RT  for  heptane  at  600  K<T<1100K 

Equilibrium  k^  K 

Reaction  AH®3oo  log  A  E  log  A'*'  E  log  A'  E" 


1  RH  +  O2  <=>  R-  H02’ 

2  R-  +  O2  <=>  R02- 

3  R02-  <=>  QOOH-  n-heptane 

63  PRF 
87PRF 
iso-octane 

4  QOOH-  +  O2  <=>  OOQOOH- 

5  OOQOOH-  =>  OQ’OOH  -t-  OH- 

6  OH-  -H  RH  =>  H2O  +  R- 

7  OQ’OOH  =>  OQ’O-  +  OH- 

8  H02-  +  HO2-  =>  HOOH  +  O2 

9  HOOH  +  M  =i.  20H-  +  M 

10  OQ’O-  =>  2  RCHO  +  RCO- 

(n-heptane) 

=>  2  RCHO  +  Rs- 

(iso-octane,  63,  and  87  PRF) 

11  QOOH-  =>C=C  +  RCHO  +  OH- 

12  RO2-  +  RCHO  =>  ROOH  +  RCO- 

13  H02-  -«-  RCHO  =i.  HOOH  +  RCO- 

14  C=C  +  HO2-  =>  Epox  +  OH- 

15  H02-  -»■  RH  <=>  R-  HOOH 

16  RO2-  -t-  RH  <=>  ROOH  +  R- 


46.4 

1.5 

46.0 

13.5 

46.0 

12.0 

0.0 

-31.0 

-1.4 

-27.4 

12.0 

0.0 

13.4 

27.4 

7.5 

0.9 

8.0 

11.9 

19.0 

11.0 

11.0 

7.5 

0.98 

11.24 

11.98 

22.4 

11.0 

11.0 

7.5 

0.59 

11.24 

11.59 

22.4 

11.0 

11.0 

7.5 

0.0 

11.24 

11.0 

22.4 

11.0 

11.0 

-31.0 

-1.9 

-27.4 

11.5 

0.0 

13.4 

27.4 

-26.6 

11.3 

17.0 

-23.5 

13.3 

3.0 

43.6 

15.6 

40.0 

-38.5 

12.3 

0.0 

51.4 

16.88 

46.0 

-17.5 

14.0 

15.0 

18.5 

14.0 

15.0 

-3.0 

14.4 

31.0 

-0.6 

11.45 

8.6 

-0.6 

11.7 

8.64 

-0.23 

10.95 

10.0 

8.0 

0.9 

8.0 

11.7 

16.0 

10.8 

8.0 

8.0 

1.1 

8.0 

11.2 

16.0 

10.1 

8.0 
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Table  2  (cont.) 

Current  Reduced  Chemical  Kinetic  Model  --  Part  2 


A.  7  Active  Species 


1. 

RCO-  2.  Rg-  3.  Rs02-  4. 

RgOOH 

5.  RgO-  6. 

RO-  7. 

ROOH 

B.  13  Reactions  (  unit:  cc.  mole.  s.  kcal  1 

Arrhenius  parameters  of  rate  constants  k  = 

Ae-E/RT 

Equilibrium 

k+ 

k- 

Reaction 

^H°30o 

log  A  E 

log  A+ 

E  + 

log  A-  E- 

17 

RCHO  +  OH-  =>  RCO-  +  H2O 

(n-heptane) 

-31.5 

1322 

0.0 

(63  PRF) 

-31.5 

13.43 

0.0 

(87  PRF) 

-31.5 

13.55 

0.0 

(iso-octane) 

-31.5 

13.57 

0.0 

18 

RCO-  M  =>  Rg-  -H  CO  -1-  M 

10.7 

16.78 

15.0 

19 

Rs•-^02<=>  Rs02- 

-31.0 

-1.4  -27.4 

12.0 

0.0 

13.4  27.4 

20 

R3O2-  =>  C=C  +  HO2- 

17.5 

11.75 

28.9 

21 

RCHO  +  Rs02-  =>  RgOOH  +  RCO-0.6 

11.53 

8.6 

22 

RH  +  Rs02-  <=>  RgOOH  +  R- 

8.0 

1.18  8.0 

1128 

16.0 

10.1  8.0 

23 

RgOOH  =>  RgO-  -1-  OH- 

43.6 

15.6 

43.0 

24 

R3O-  -1-  O2  =>  Rs'O  -•-  H02- 

-26.5 

10.6 

2.14 

25 

C=C  -1-  OH-  =>  2  OXY  -1-  OH 

-75.5 

12.72 

-1.04 

26 

ROOH  =>  RO-  -t-  OH- 

43.6 

15.6 

43.0 

27 

RO-  =>  Rg-  -h  RCHO 

-10.0 

13.3 

15.0 

28 

R02-  =>  C=C  +  HO2- 

4.0 

9.85 

23.0 

29 

RO2-  =>  ether  -1-  OH-  (n-heptane) 

-25.0 

9.48 

18.0 

(63  PRF) 

-25.0 

928 

18.0 

(87  PRF) 

-25.0 

9.18 

18.0 

(iso-Octane) 

-25.0 

8.78 

18.0 

16 


